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Abstract 

The dimensionally regularized massless on-shell double box Feynman 
diagram with powers of propagators equal to one is analytically evalu- 
ated for general values of the Mandelstam variables s and t. An explicit 
result is expressed either in terms of polylogarithms Li a (—t/s), up to 
a = 4, and generalized polylogarithms S at b(—t/s), with a = 1,2 and 
b = 2, or in terms of these functions depending on the inverse ratio, 
s/t. 
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1 Introduction 



The massless double box diagram shown in Fig. 1 enters many important physical 
observables, e.g., amplitudes of the Bhabba scattering at high energies. An experience 
shows that master diagrams, i.e. with all powers of propagators equal to one, are most 
complicated for evaluation. In the massless off-shell case, the master double box 
Feynman integral has been analytically evaluated in strictly in four dimensions. 
It is the purpose of the present paper to evaluate it analytically on shell, i.e. for 
p 2 = 0, i = 1,2,3,4, in the framework of dimensional regularization JJ], with the 
space-time dimension d = 4 — 2e as a regularization parameter. 




To do this, we start, in the next Section, from the alpha-representation of the 
double box and, after expanding some of the involved functions in Mellin-Barnes 
(MB) integrals, arrive at a five-fold MB integral representation with gamma func- 
tions in the integrand. We use, in Sec. 3, a standard procedure of taking residues 
and shifting contours to resolve the structure of singularities in the parameter of di- 
mensional regularization, e. This procedure leads to an appearance of multiple terms 
where Laurent expansion in e becomes possible. The resulting integrals in all the MB 
parameters but one are evaluated explicitly in gamma functions and their derivatives. 
In Sec. 4, the last MB integral is evaluated by closing an initial integration contour 
in the complex plane to the right, with an explicit summation of the corresponding 
series. A final result is expressed in terms of polylogarithms Li a (— t/s), up to a = 4, 
and generalized polylogarithms S a ^{—t/s), with a = 1,2 and b = 2. Starting from 
the same one- fold MB integral and closing the contour of integration to the left, we 
obtain a similar result written through the same class of functions depending on the 
inverse ratio, s/t. Furthermore, we obtain, as a by-product, an explicit result for the 
backward scattering value, i.e. at t = —s, of the double box diagram. 

2 From momentum space to MB representation 

The massless on-shell double box Feynman integral can be written as 




d d kd d l 

(k 2 + 2 Pl k)(k 2 - 2p 2 k)k 2 {k - l) 2 (l 2 + 2 Pl l){l 2 - 2p 2 l){l - pi -p 3 ) 2 



1 
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K(t/s,e) 



where s = (pi +P2) 2 , t = (pi +P3) 2 , and k and I are respectively loop momenta of 
the left and the right box. Usual prescriptions, k 2 = k 2 + iO, —s = —s — iO, etc are 
implied. We have pulled out not only standard factors that arise when integrating in 
the loop momenta but also a factor that makes the resulting function K depend on 
the dimensionless variable, x = t/s. 

The alpha representation of the double box is straightforwardly obtained: 
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where 



D = (ati + a 2 + a 7 )(a 3 + a 4 + a 5 ) + a 6 («i + a 2 + a 3 + a 4 + a 5 + a 7 ) , (3) 
A = aia 2 («3 + "4 + oc 5 ) + a 3 a 4 (ai + ct 2 + 017) + a 6 (a;i + a 3 ){a 2 + a 4) . (4) 

As it is well-known, one can choose a sum of an arbitrary subset of a.i , i — 1, . . . , 7 in 
the argument of the delta function in @. We choose it as 5 a i ~ l) an d change 

variables by turning from alpha to Feynman parameters 

"3 = "356, "5 = "35(1 - 6), "1 = "176, "7 = "17(1 - £3), 

«35 = 662, a 4 = £5(1 - 62), «i 7 = (1 - £ 5 )£ 4 , a 2 = (1 - 6)(1 - 6)- (5) 
to obtain the following parametric integral: 

K(x, e) = -r(3 + 2e) / da 6 / d& . . . / d£ 5 66f(l - 2 

JO JO JO 

x(a 6 + e 5 (l-e 5 )) 1+3e g- 3 - 2 % (6) 

where 

Q = za 6 (l -6)6(1-6)6(1- 6)6 

+6(1 - 6) [666(1 - 6) + (1 - 6)66(1 - &)] 

+« 6 [666 + (1 - 6)66][6(i - 6) + (1 - 6)(1 - 6)] ■ (7) 

We are now going to apply five times the MB representation 

I 11 r+ioo y w 

dw—— T(u + w)T(-w) , (8) 



(X + Y) v Y{u) 2m J-ioc X"+ 

where the contour of integration is chosen in the standard way: the poles with the 
T(. . . + u>)-dependence (let us call them infrared (IR) poles) are to the left of the 
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contour and the poles with the T(. . . — u>)-dependence (ultraviolet (UV) poles) are to 
the right of it. 

First, we introduce a MB integration, in w, using decomposition of the function Q 
with Y as the first line in (0). We introduce a second MB integral choosing as X the 
term with in the rest part of Q. After that we can take an integral in a§ in gamma 
functions. The next three MB integrations, in Z\,z 2 and z 3 , are to separate terms in 
the following three combinations: [£5616 + (1 - 65)6361], [Cs(l - 62) + (1 - 6s)(l - 64)] 
and [65616(1 - 62) + (1 - 65)6364(1 - 64)]. 

All the integrals in Feynman parameters are then taken explicitly in gamma func- 
tions. Finally, we perform the change of variables Z2 = w 2 — Z\ — 1, z 3 = w 3 — z\ — 1 
and arrive at the following nice 5- fold MB integral: 

1 1 f 

K(x,e) = ; — / dwdw 2 dw 3 dzdz 1 x w+1 

v ' ; r(-l-3e) (2ttz) 5 J 231 

xT(l + w) 2 T(-w)T{w 2 )T(-l -2e-w- w 2 )T(w 3 )T(-l - 2e - w - w 3 ) 
T(l - w 2 + zi)r(l - w 3 + z l )T(e + w + w 2 + w 3 - Zi)T(—Zi) 
T(l + w + w 2 + w 3 )T(-l — 4e — w — w 2 — w 3 ) 
xT(l - e + z)r(2 + 2e + w + w 2 + z- z 1 )T(2 + 2e + w + w 3 + z - z x ) 

r(-2 - 3e - w - w 2 - w 3 + zi - z)Y{z 1 - z) 
X r(3 + 2e + w + ^) ' U 

One can interchange the order of integration in an arbitrary way. For each order, 
the rules of dealing with poles are as formulated above. Note that if we have a 
product T(a + v)Y{b — v), for some integration variable v = w, w 2 , w 3 , z, Z\ with a and 
b dependent on other variables, then the integration in v produces a singularity of 
the type T(a + b). 



3 Resolving singularities in e 

Since it looks hopeless to evaluate our MB integral for general e let us try to obtain 
a result in expansion in e up to the finite part. There is a factor 1/T(— 1 — 3e) 
proportional to e when e tends to zero. Representation (|9|) is therefore effectively 
4-fold because to generate a contribution that does not vanish at e = we need to 
take a residue at least in one of the integration variables. None of the integrations can 
however immediately produce an explicit e-pole. Let us first distinguish the following 
two gamma functions 

T(e + w + w 2 + w 3 - zi) r(— 2 — 3e — w — w 2 — w 3 + z x — z) 

that are essential for the appearance of the poles. 

We can write down the integral in Z\ as minus residue at the point z\ = e + w + 
w 2 + w 3 (where the gamma function r(e + w + w 2 + w 3 — Zi) has its first pole which 
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is UV, with respect to z\) plus an integral with the same integrand where this pole 
is IR. We can similarly write down the integral in z as minus residue at the point 
z = —2—3e—w—w 2 —w 3 +Zi (where the gamma function T(— 2— 3e—w—w 2 —w 3 +zi — z) 
has its first pole which is UV, with respect to z) plus an integral with the same 
integrand where this pole is IR. 

As a result we decompose integral (|) as K = K 00 + K m + K w + K n where Ku 
corresponds to the two residues, K w to the residue in z and the integral in Z\ with 
the opposite nature of the first pole of T(e + w + w 2 + w 3 — zi), etc. For example, the 
contribution Kn is given by the following 3-fold integral: 



K n (x,e) = / dwdw 2 dw 3 x w+1 T(l + w)T( 



w 



(2vri 

xr(w 2 )r(-e - tw 2 )r(l + e + w + tw 2 )r(-l - 2e - w - w 2 
xr(w 3 )r(-e - w 3 )T(l + e + w + w 3 )T(-l - 2e - w - w 3 
T(2 + 3e + w + w 2 + w 3 )T(-e - w - w 2 - w 3 ) 



(10) 



T(l + w + w 2 + w 3 )T(— 1 — 4e — w — w 2 — w 3 ) 

This contribution is in turn decomposed, in a similar way, as K u = J2i,j=o,i,2^mj- 
Here the value i — 1 of the first index denotes the residue in w 2 at the point w 2 = 0. 
The value i = 2 denotes the residue in w 2 at w 2 = — 1 — e — w of the integrand where 
the first pole of T(w 2 ) is UV rather than IR. Finally, % = means that both above 
poles are IR. The second index similarly refers to the integral in w 3 . 
In particular, we have 

2 1 I A„„„W+1 



T(l + e + w) 2 T(2 + 3e + w)T(-l - 2e - w) 2 T(-e - w)V(-w) 
X r(-l-4e-«;) ' ( } 

Knn(x,e) = K im (x,e) = F(1 + 2c f ( " c) -L f dwx ^ 

1 (— 3eJ 2m J 

r(l + w) 2 T(l + t + w) 2 T{-\ - 2e - w)T(-e - w)T(-w) 

F(2 + e + w) 



Kimix, e) = _r ( _e ) 2 ^~~ / du7X 
r(l + w) 3 T(l + e + w) 2 T{-e - w) 2 T{e - w)T(-w) 
X T(2 + e + w)T(l-2e + w)T(-l-2e-w) ' ^ J 

The next contribution is 

Kio(x,e) = - r( _ 1 1 _ — ^ J dwdw 2 dw 3 d Zl x w+1 T{l + w) 2 T( 

xT(w 2 )T(-e - m 2 )r(-l -2e-w- ^ 2 )r(w 3 )r(-e - w 3 )T(-l - 2e - w - w 3 ) 

r(2 + 3e + w + w 2 + m 3 )r(l - w 2 + ^)r(l - w 3 + zi) 



-W 



X- 



T(l + w + w 2 + w 3 )T(— 1 — 4e — u; — w 2 — w 3 ) 
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x r(-l - 4e - w - w 2 - w 3 + z 1 )T(e + w + w 2 + w 3 - gi)T(-Zi) 

r(l — e — w 2 — w 3 + Zx) 

where the first pole of T(e + w + w 2 + w 3 — Z\) is IR, with respect to z±, rather than 
UV. We further decompose this contribution by changing the nature of the first pole 
of T( — 1 — 4e — w — w 2 — w 3 + zi) in z±. We obtain K w = K wo + K W i, where the new 
index 1 corresponds to the residue and has the form 

ivioifx, e) = - / dwdw 2 dw 3 x + 



(2mf J r(2 + 3e + w) 

xr(w 2 )r(-e - w 2 )Y(2 + 4e + w + iy 2 )r(-l -2e-w- w 2 )T(w 3 ) 
T(-e - w 3 )T{2 + 4e + w + w 3 )T{-l - 2e - w - w 3 )T{2 + 3e + w + w 2 + w 3 ) 
X T{l + w + w 2 + w 3 ) 

Each of the contributions K wo and i^ioi is then decomposed using the change of the 
nature of poles w 2 = and w 3 = 0. We obtain K 10 j = i^iojoo+-^ioioi+-^iOjio+-^ioiii) 
for j = 0,1. Here the value i = 1 of the last index denotes the residue in w 3 at the 
point w 3 = and the i = an integral where the first pole of r(u> 3 ) is considered UV. 
The second index from the end similarly refers to T(w 2 ). For example, 

#10111 (a, e) = -r(-e) 2 -^ / dwx w+1 T(2 + 4e + w) 2 T(l+w)T(-l-2e-w) 2 T(-w). 

2m J 

(16) 

Then we have 

K 01 (x,e) = - r( _ 1 1 _ 3e) ( ^. )4 J dwdw 2 dw 3 dzx w+1 T(l + w) 2 T(-w) 

xT(w 2 )T(l + e + w + w 2 )T(-l - 2e-w- w 2 )T(w 3 )T(l + e + w + w 3 ) 
T(— 1 — 2e — w — w 3 )T(— e — w — w 2 — w 3 )T(l — e + z) 



T(l + w + w 2 + w 3 )T(—l — 4e — w — w 2 — w 3 ) 
T(2 + e-w 2 + z)T(2 + e-w 3 + z)T(e + w + w 2 + w 3 - z)T(-2 -2e-z) 



T(3 + 2e + w + z) 



(17) 



where the first pole of T(— 2 — 2e — z) is IR, with respect to z, rather than UV. Using 
the change of variables w 2 — > — 1 — 2e — w — w 2 , w 3 — >• — 1 — 2e — w — u>3, — > 
— e — w — W2 — w 3 + z in K w we see that Xoi = Kio- Finally, the contribution K Q0 is 
similarly decomposed: K 00 = K 0000 + K 0001 + K oow + # on- 

Now we observe that, in each of the obtained contributions, the only additional 
(with respect to explicit gamma functions depending on e) source of the poles in 
e is the last integration, in w, where the first (UV) pole of the gamma function 
T(— 1 — 2e — w) glues with an IR pole of T(l + w) or T(l + e + w) when e — > - 
see such examples in (|Il"|--|i~3"|) and flUf). Therefore we further decompose each of the 



contributions into two pieces: minus residue at the point w = — 1 — 2e plus an integral 
where we can integrate in the region — 1 < Kew < 0. In each of these pieces, we now 
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can expand an integrand in a Laurent series in e up to the finite part. In particular, 
no poles in e arise in -Koooo so that it it zero at e = because of the overall factor 
l/r(-l-3e). 

We collect separately the pieces from these last residues and from the last inte- 
gration at — 1 < Rew < 0. The first collection gives the leading order term in the 
expansion of the double box in the limit t/s — >• while the second collection involves 
the rest of the terms of this expansion. A remarkable fact is that, in all these multiple 
contributions, the integrations in w 2 , tu 3 , z, Z\ can be performed analytically, with the 
help of the first and the second Barnes lemmas 



I r+ioo 

— / dw r(Ai + w)T(X 2 + w)T(X 3 - w)T(X 4 - w) 

2m J-ioo 

r(Ax + A 3 )r(Ax + A 4 )r(A 2 + A 3 )r(A 2 + a 4 ) 

T(Ai + A 2 + A 3 + A 4 ) 

r(Ax + w)r(A 2 + w)r(A 3 + w)r(A 4 - w)r(A 5 - w) 

T(Ai + A 2 + A 3 + A 4 + A 5 + u>) 
= T(X 1 + A 4 )r(A 2 + A 4 )r(A 3 + A 4 )r(Ax + A 5 )r(A 2 + A 5 )r(A 3 + As) 
r(Ai + A 2 + A 4 + A 5 )r(Ai + A 3 + A 4 + A 5 )r(A 2 + A 3 + A 4 + A 5 ) 
and their corollaries. These are two typical examples of such corollaries: 
1 f+ io ° ^ r(Ai + w)T(X 2 + w) 2 T(-X 2 - w)T(X 3 - w) 



(18) 



1 /•+< 
2m J-ic 



(19) 



2-ni J-ioo T(Ai + A 2 + A 3 + w) 

r(A x - A 2 )r(A 2 + A 3 ) W (Ai + A 3 ) - (A 2 + A 3 )] 



(20) 



r(Ax + A 3 ) 

where the pole w = — A 2 is considered IR while other poles are treated in the standard 
way, and 

— / dw T(l + w)T(w)r(-w)T(-l - w)1>(l + w) 2 = ^- + 6 7E C(3) + - . 

2m J-i/2-ioo 3 45 

(21) 

Here 7e is the Euler constant, ip(z) the logarithmical derivative of the gamma func- 
tion, and ((z) the Riemann zeta function. 

After taking these integrations and summing up the resulting contributions into 
the two above collections we obtain the following result 

K(x,e) = K ot (x, e) + K lt (x, e) + o(e) , (22) 

„ , s 4 5 In a; / , 9 5 ,\ 1 

K ot (x,e) = -- + —^-[2ln 2 x--n 2 j- 

_ ( ^ ln 3 x + 11^2 ln _ 65 C(3) A I 1 4 2 2 _ 88 29 4 

V3 2 3 v V e 3 3 w 30 v ' 

X lt x,e = - / — —T(l + w) 3 T(-w) 3 
m J 1 + w 

(24) 



+ w 
I 5 

+ 3^(1 + w) — 4^(— iu) — 7e 



1 + iu 
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Let us stop for a moment and observe that this result provides, in a very easy 
way, not only numerical evaluation of the double box diagram for general values of s 
and t but also asymptotic expansions in the limits t/s — > and s/t — > which are 
obtained by taking series of residues respectively to the right or to the left. 



4 Evaluating the last MB integral 



The last MB integration, in (p4|) , is performed analytically by taking the sum of the 
residues at the points w = 0, 1,2, . . . and summing up the resulting series. In this 
last step, we use, in particular, summation formulae derived in ||. Here is the final 
result: 



K lt (x,e 



2Li 3 (-x) - 2 lnxLi 2 (-x) - (in 2 x + vr 2 ) ln(l + 



x : 



-4 (S 2 ,2(-x) - \nxS h2 {-x)) + 44Li 4 {-x) - 
+2 (ln 2 x + 21nxln(l + x) H 7r 2 ) Li 2 (-x 



4 (ln(l +x) + 61nx) Li 2 



-X) 



+ (in 2 x + vr 2 ) ln 2 (l + x) - - (4 In 3 x + 5vr 2 In x - 6C(3)) ln(l + 

where Li a (z) is the polylogarithm [Q] and 

(_l)o+6-i i \ n a -\t)\n\l- zt) 



X 



(25) 



S a ,b(z) 



(a-l)\b\ 



-dt 



(26) 



a generalized polylogarithm introduced in [Q. Note that any (generalized) polylog- 
arithms involved can be expanded in a Taylor series at x = with the radius of 
convergence equal to one. 

We can similarly close the integration contour to the left and obtain a result in a 
form of functions depending on the inverse ratio, y — 1/x: 



K(x,e) 
K 0s (l/y,e) 



K 0s (x, e) + K ls (x, e) + o(e) 



4 5 \ny 



65 



e 

1 1 



21n 2 y-^)- 
76 



S3 



+ -vr 2 \ny + — C(3) - + -vr 2 In 2 y + —((3) lny - — vr 



90 



(27) 



(2? 



2Li 3 (-y) - 2 lnyLi 2 (-y) - (in 2 y + tt 2 ) ln(l + y) 



-4 OM-j/) - lny^C-y)) - 36Li 4 (-y) - 4 (ln(l + y) - 5 Iny) Li 3 (-y) 
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-2 [\n 2 y-2\ny ln(l + y) + — tt 2 J Li 2 (-y) 
+ (in 2 y + vr 2 ) ln 2 (l + y) + 2 (in 3 y + |vr 2 In y + 2((3)) ln(l + y) . 



(29) 
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As a by-product, we obtain an explicit result for the backward scattering value of 
(|l|), i.e. at t — —s, 



'_ s )3+2e 



4 9vr 2 53C(3) 22tt 4 / 5 25vr 2 148C(3V 

— r r 1 7T« -r 

e 4 2e 2 3e 9 1 e 3 6e 3 



(30) 



The presented algorithm is applicable to massless on-shell box Feynman integrals 
with any integer powers of propagators. 
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